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1.  INTRODUCTION 


Energy  exchange  and  upscattering  in  ion-ion  collisions  as  they  pass  through  the 
core  of  a  Polywell*"^/HEPS  device  have  been  analyzed  previously  by  Lovberg*  and 
Bussard^.  One  result  of  these  studies  is  a  simple  equation  describing  the  energy  dispersion 
introduced  over  a  time  t^  by  such  collisions.  This  is  given  by^ 

<Ej(t^)>  =  x(2nvt/\Ze)^  (1) 

Here  <E^(t^)>  is  the  total  energy  spread  introduced  in  an  ion  population  by  small  angle 

collisions  in  the  core,  after  repeated  passes  through  the  core  by  the  subject  ions.  This 

*  • 

spread  both  increases  and  decreases  the  energy  of  ions  from  their  initially— monoenergetic 
level  at  the  well  center.  Taking  fE  as  the  mean  up  OR  down  dispersion  of  energy  from 
in— core  scattering  collisions,  gives  26E  =  ±  <Ej(t^)>.  Defining  a  function  f  =  (6v/v)  as 
the  fractional  up  or  down  spreading  in  velocity  space  gives  {6E/E)  =  2f,  so  that  ^E  =  2fE. 
With  this,  eq.  (1)  becomes 

4fE  =  7r(2nvt^)®-^(Ze)^  (2) 


2.  UPSCATTERING  TIME 

This  can  be  solved  for  t^  as  the  time  required  to  3rield  an  energy  upscatter  of  2fE  or 
a  fractional  velocity  upscatter  of  f.  In  order  for  this  to  be  conect,  the  product  (nv)  must 
be  the  integrated  average 

<nv>  =  (l/R)jo“(rMr)dr 
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over  the  complete  transit  path  of  the  ion.  This  is  a  measure  of  the  scattering  center 
frontal— area-loading-<iensity  seen  by  the  particles  as  they  circulate  across  the  device. 

For  purposes  of  assessing  energy  broadening  (upscattering  and  downscattering)  of 
the  convergent  ion  flow  the  scatterings  that  are  effective  are  only  those  in  which  the 
scattered  ions  have  significant  energy  themselves.  Scatterings  that  occur  among  low  energy 
ions  at  the  periphery  of  the  system  (or  at  or  near  the  ion  injection  point  within  the  system) 
will  not  lead  to  energy  broadening  (but  may  have  other  effects;  these  are  treated  in  another 

n  * 

EMC  technical  note)  . 

Thus  it  is  sufficient,  in  estimating  <nv>  to  take  the  ion  speed  as  approximately 
constant  at  a  mean  core  speed  v^  in  the  (extensive)  region  of  the  bottom  of  the  potential 
well,  about  the  system  center.  This  allows  the  product  <nv>  to  be  written  as 
<nv>  ss  <n>v^.  Over  this  central  region,  the  average  ion  density  is  comprised  of  two 
terms,  that  due  to  core  ions  within  r  <  r  and  that  due  to  ions  outside  the  core 
convergence-limited  radius  r  =  r^. 

In  the  core  region  the  ion  density  is  assumed  constant  at  n^within  the  radius  r^. 
Outside  r^  the  results  of  simple  analyses  and  extensive  KXL/EKXL  computations  show 
that  the  density  varies  as  the  inverse  square  of  the  radial  distance  from  the  core,  out  to  a 
radius  r  =  r^  at  which  the  potential  begins  to  rise  towards  the  outer  edge  of  the  system. 
This  point  is  found  to  be  generally  in  the  range  0.3  <  <rj>  <  0.7,  and  thus  approaches 
that  of  the  ion  injection  origin.  The  density  variation  of  importance  here  can  then  be 


Outer  region  ion— ion  collision!  may  result  in  increased  losses  by  deflection  scattering,  and  reduced 
transverse  momentum  by  isotropisation  of  outer  peripheral  low— energy  ions. 
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written  simply  as  n(r)  =  n^(r^/r)  .  Using  these  the  path-averaged  ion  density  becomes 


<n>  =  2(r^/R)n^  (3) 

With  this,  and  writing  =  (v^/v^)(2E/m.)®‘^,  where  is  the  ion  speed  at  the 
core  at  energy  E,  eq.  (2)  can  be  solved  for  t^  as 


t 

c 


4f^E^-®  fR 
7r^(Ze/n^  r^ 


(4) 


This  equation  gives  the  time  for  ffi/E  =  2f  broadening  of  ion  energy  for  collisions  at 
and  around  an  energy  E.  But  the  ion  losses  of  concern  are  due  to  upscattering  (rather 
than  downscattering),  thus  the  effect  of  reduced  collision  cross-sections  in  higher  energy 
collisions  must  also  be  taken  into  account  in  its  estimation.  As  shown  by  Rosenbluth, 
et  aH’5,  and  noted  previously*,  the  time  required  to  upscatter  to  the  maximum  velocity 
v^  =  v^(l+f)  is  increased  by  a  factor  of  about  one-half  of  the  cube  of  the  ratio  of  this 
velocity  to  the  initial  (mean)  ion  velocity  v^;  (l/2^(vj/v^)^.  Since  this  is  true  for 
upscattering  to  all  velocities  between  these  two,  the  upscattering  time  factor  averaged  over 
all  of  the  particles  is  just  this  cubic  dependence  weighted  by  the  Maxwellian  distribution, 
integrated  over  v^  <  v  <  v^  This  gives  a  weighted  average  factor  of  approximately 
(l/2)(l+f)  ,  which  must  be  included  in  the  time  estimate  of  eq.  (4).  This  gives 


t 

C 


ir^  (Ze)^n^<r^> 


(5) 


where  <r^>  =  (r^/R)  is  the  convergence  ratio  of  the  system 7,  and  the  ion  mass  is  taken  as 
m.  =  Am^  where  A  is  ion  mass  number  and  m^  is  the  proton  mass.  Taking  units  of  eV  for 


4 


£  and  all  others  in  cgs  gives 


t 

c 


0.71E7 


0.5 


n  <r  >  Z 

c  c 


(6) 


For  the  example  case  given  earlier  by  Krall  and  Rosenberg*,  with  Z  =  1,  A  =  2, 
<r^>  =  2E-2,  v^/v^  »  1,  f  =  0.1,  E  =  1E4  eV,  and  n^  =  lE12/cm^,  this  equation  gives 
t^  =  6.1  sec.  Two  other  cases  they  examined  later  were  for  <r^>  =  lE-2,  and  f  =  0.5, 
with  other  parameters  as  above,  and  for  this  same  case  but  with  higher  density, 
n^  =  lE18/cm^,  and  energy,  E  =  1E5  eV.  For  these  eq.  (6)  gives  the  upscattering  times  as 
t^  =  565  sec  and  t^  =  17.9E— 3  sec,  respectively.  All  of  these  calculated  values  agree  closely 
with  results  of  their  Fokker— Planck  analysis  for  these  examples.  This  simple  formula  thus 
can  be  used  to  model  collision  times  for  upscattering. 


3.  COLLISIONAL  LOSSES 

Such  core  collisionally-driven  upscatter  of  ions  will  lead  to  expansion  of  their 
orbital  motion  to  regions  of  higher  magnetic  field  beyond  that  (outside)  of  their  injection 
point.  If  the  energy  upscatter  is  great  enough,  the  ions  can  leave  the  well  entirely.  This  is 
the  case  considered  here. 


If  upscatter  is  to  lesser  energies,  circulation  to  larger  radii  may  result  in  greater 
transverse  energy /momentum  generation  by  gyro  motion  within  the  surface  field  regions 
which,  in  turn,  can  result  in  an  increase  in  the  ion  convergence  ratio  at  the  core,  and  thus 
yield  lower  core  densities  and  poorer  performance.  This  latter  case  is  considered  in  another 
EMC  technical  note. 

For  this  case  of  upscatter  to  direct  escape  energies,  the  collision  time  t^  (eqs.  5,6) 


5 


can  be  used  to  formulate  a  loss  term  in  the  algorithm  for  ion  edge  density  (n^.  at  r  =  R) 
time-dependence  used  in  the  EKXL  code.  Using  the  forms  derived  by  King  and  Bussard®, 
and  carrying  out  the  algebra,  gives 


Here,  as  before®,  average  ion  density 

O 

is  <n.>  =  N./(47r/3)R  ,  and  is  the  anode  height  control  function  that  limits  the 
time-incremental  density  change  used  in  the  EKXL  computation,  and  The 

term  is  the  ion  current  recirculation  ratio  in  the  system,  as  set  by  ion  loss  mechanisms 
other  than  direct  escape  by  collisional  upscattering  (e.g.  by  particle  cusp  losses). 
Implementation  of  this  algorithm  has  been  made  in  the  EKXL  code,  v.2.0.  Calculational 
runs  made  with  this  code  show  that  ion-ion  core  upscattering  collisions  have  virtually  no 
effect  in  systems  of  interest  at  central  densities  below  about  lE18/cm^.  Since  this  density 
regime  is  already  high  enough  to  give  interesting  fusion  output  and  system  gain  relative  to 
electric  drive  power,  the  phenomenological  effect  of  this  model  of  upscattering  is  not 
dominant  for  most  systems  of  practical  interest. 
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